A molecular dynamics simulation of polymer crystallization 
from oriented amorphous state 
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Molecular process of crystallization from an oriented amorphous state was reproduced by molec- 
ular dynamics simulation for a realistic polyethylene model. Initial oriented amorphous state was 
obtained by uniaxial drawing an isotropic glassy state at 100 K. By the temperature jump from 100 
K to 330 K, there occurred the crystallization into the fiber structure, during the process of which 
we observed the developments of various order parameters. The real space image and its Fourier 
transform revealed that a hexagonally ordered domain was initially formed, and then highly ordered 
crystalline state with stacked lamellae developed after further adjustment of the relative heights of 
the chains along their axes. 

PACS numbers: 61.50.-f, 61.41.+e, 83.10.Rs 



The molecular dynamics (MD) simulation has recently 
come to be recognized as a very promising tool to un- 
derstand crystallization in systems of chain molecules, 
since they can directly provide us with molecular trajec- 
tories that are experimentally hard to access. In an early 
paper Rigby and Roe successfully reproduced a highly 
ordered lamella structure from a bulk amorphous state 
of 71-decane The work was followed by studies that 
elaborated on detailed crystallization process in bulk n- 
alkanes of various length ^. Crystallization in thin 
films of ?i-alkanes was also investigated using MD tech- 
nique 1^, [sj. Recent advances have come to enable the 
quantitative estimation of the growth rate vs. crystal- 
lization temperature |^ . 

Chain folded crystallization of polymers has also been 
investigated by molecular simulations. Kavassalis and 
Sundararajan first reported the MD studies on the col- 
lapse and folding of a single polyethylene molecule in vac- 
uum 0, while Muthukumar and coworkers investigated 
a primary nucleation and subsequent crystal growth in 
solution by Langevin dynamics simulation ^. The 
molecular process of secondary nucleation and growth of 
lamella crystals has also been a subject of great interest 
and several investigations have been reported [p|, p^, pd]|. 
We can notice in very recent literatures a surge of interest 
in reproducing chain folded crystallization from the melt 
and solution both during the primary nucleation |l^] 
and the growth |]l3| , though at present we can only treat 
relatively short polymer chains of about hundreds of seg- 
ments. 

Crystallization of much longer chains in the bulk sys- 
tem is expected to require extraordinary long simulation 
time even by use of present day supercomputers, espe- 
cially when we adopt realistic molecular models. We have 
already reported an MD simulation of polymer ordering 
in the bulk polyethylene We could there observe 

only the local ordering in chain packing, and we discussed 



the results within the context of recent experimental and 
theoretical investigations on precursor phenomena prior 
to crystalhzation |l5[ [6[ 

The polymeric material needs very long time until the 
onset of crystallization and can be easily quenched into 
amorphous glassy state with isotropic chain orientation. 
By drawing the amorphous sample below the glass tran- 
sition temperature, we obtain the oriented amorphous 
state in which the chains run parallel to the draw di- 
rection on the average. We can then expect to accel- 
erate the crystallization if we use such oriented amor- 
phous sample, instead of isotropic melt JT^ , as an ini- 
tial state of MD simulation. In fact, recent experiments 
revealed that the crystallization rates of polymers de- 
pend strongly on the degree of the chain orientation [ p^ . 
Moreover, it will be found that the oriented amorphous 
state is very convenient in investigating the spontaneous 
development of fiber structures with stacked lamellae or 
the detailed molecular level structure of the amorphous 
phase between lamellae . 

In this Letter we report results of our recent simulation 
of polymer crystallization from bulk oriented amorphous 
state. The molecular model we adopted was essentially 
the same as that of our previous work . We considered 
the energies of bond stretching, bond angle bending, and 
torsion of the united atom model of polyethylene (PE). 
The truncated Lennard-Jones 12-6 potential was used to 
reproduce van der Waals interactions between the atoms, 
where a rather short cutoff radius of 0.85 nm was used 
to improve computational efficiency. Bulk polymer sys- 
tem was made with a single linear PE composed of 5000 
united atoms being placed in the MD cell under periodic 
boundary condition. Newton's equations of motion were 
integrated by the leap- flog algorithm with a time step of 4 
fs, and the temperature was controlled every 10 steps by 
the ad-hoc velocity rescaling, while the pressure was ad- 
justed every 100 steps by the loose-couphng method pl| . 



2 




FIG. 1: Schematic picture of the preparation of oriented 
amorphous sample. 

The parameter M used in the pressure control was set to 
be 13.26 Pa-s/ni, and the shape of the MD cell was kept 
orthorhombic. 

First we prepared the isotropic melt as follows. An ini- 
tial configuration of the chain, with fixed bond lengths 
and bond angles, was generated using the technique de- 
veloped by Theodorou and Suter and then a short 
MD simulation of 10 ps under constant volume and tem- 
perature (600 K) was made to reduce the overlap of 
atoms. This initial state was then fully relaxed for 1 ns 
at constant temperature (600 K) and pressure (1 atm). 

Next we generated the oriented amorphous state by 
drawing the isotropic amorphous sample well below the 
glass transition temperature. The glassy amorphous 
sample was drawn at lOOK along the x-axis at a con- 
stant rate of 500 % /ns up to a final draw ratio of 400 % 
(Fig.l), where the method of uniaxial drawing was the 
same as that used by Ogura and Yamamoto p3| . During 
drawing, the temperature was kept constant at 100 K 
and the lateral widths of the cell were allowed to shrink 
along the directions perpendicular to the draw axis un- 
der constant pressure (1 atm). After attaining the final 
draw ratio of 400%, further relaxation was carried out for 
1.5 ns under constant pressure (1 atm) and temperature 
(100 K). 

Crystallization of the oriented amorphous state was 
then simulated for 30 ns, under constant temperature 
(330 K) and pressure (1 atm). Figure 2 shows the changes 
with time of various quantities of the system. The inter- 
nal energy slowly decreases until 10 ns (Fig. 2 (a)), but 
shows the steep reduction around 12 ns suggesting an 
onset of ordering. The specific volume also decreases al- 
most lineally up to around 14 ns (Fig. 2 (b)), where 
the MD-cell shows clear initial shrink within 4 ns and 
then expands along the draw axis (Fig. 2 (c)) at the 
most pronounced rate around 12 ns. Around 15 ns, the 
MD-cell dimensions nearly come to saturate. It might 
seem strange that the specific volume increases again af- 



ter 15 ns, but later this will be found an artifact due to 
the imposed periodic boundary condition and very slow 
relaxation in the well developed ordered state. 

PE molecule has the lowest energy conformation trans 
and the second lowest gauche^. Since the molecule has 
all trans conformation in the crystal, we evaluated the 
chain ordering by use of the trans fraction and the mean 
length of the trans sequence; the trans and the gauche 
bonds are here defined as the bonds whose torsion angles 
distribute around 0° and ±120°, respectively, within the 
range of 60°, whereas the length of the trans sequence 
is expressed by the number of the trans bonds in the se- 
quence. The transition rates in Fig. 2(d) slowly decrease 
until around 10 ns, but then they abruptly fall around 12 
ns until the transitions nearly halt around 17 ns. Both 
the trans fraction and the mean length of the trans se- 
quences in Figs. 2 (e) and 2 (f) increase steadily and 
show a marked jump up around 12 ns until they satu- 
rate completely around 17 ns. The final state obtained 
after about 20 ns is thus considered to be made predom- 
inantly of the trans segments and highly ordered with 
very infrequent internal rotations between the trans and 
gauche. 

In order to investigate the development of orientational 
order of the molecule, we examined the orientation of 
the C-C-C segments and the correlation volume that is 
defined as a measure of parallel alignment of the segments 
||l|. The orientations of the C-C-C segments along the 
three axes are evaluated by the mean squared cosines 
< cos^ 6x >, < cos^ 9y >, and < cos^ 6z >, where 9s 
are the angles between the orientation vector spanning 
every other carbons of the C-C-C segment and the axes. 
The correlation volume is defined as 

< Vc >=< / 4:'Kr'^g{r)s{r)dr >, 
Jo 

where g{r) is the radial distribution function of the cen- 
ter of mass for the segment, and s(r) is the orienta- 
tion correlation function between two segments defined 
as s(r) = 1/2(3 < cos^ 9ij{r) >ij), where 9ij{r) is the 
angle between the z-th and the j-th orientation vectors 
separated by r, and the average is taken over all pairs 
within the shell r ^ r + Ar, and the upper limit of the 
integration rmax — 1-6 nm. The chain orientation along 
the draw axis and the correlation volume were initially 
unstable, but started to increase around 4 ns until they 
saturate around 20 ns (Figs. 2(g), 2(h)). The rapid in- 
crease in the correlation volume around 12 ns indicates 
the spontaneous orientation and the parallel alignment 
of the segments. 

The structural changes described above were analyzed 
by use of real-space image and its Fourier transform, the 
structure function. Three dimensional structure function 
is defined by 

S{(1) ^ J J rfrc?r'<p(r)p(r')>exp(-iq-(r-r')), 

where q is the wave vector, r and r' are the position 
vectors, N is the total number of atoms. The atomic 
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FIG. 2: Changes in the various quantities with time: (a) the internal energy, (b) the specific volume, (c) the side length of 
the MD cell, (d) the conformational transition rate, (e) the trans fraction, (f) the mean length of the trans sequences, (g) the 
degree of the orientation, (h) the correlation volume. 





FIG. 3: Snapshots (left), and the structure functions (right) averaged around the draw axis, S2o{q\\,q±). 
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density at the position r, p{r), is defined by the relation 
p{r) = X^jLi '^(i' ~ Tj), where Vj represents the position 
of the j-th atom. Here we consider the structure func- 
tion averaged around the draw axis, 520(911 , where 
q\ I and q± are the components of the wave vector along 
the draw axis and perpendicular to it, respectively. The 
intensity profile at 10 ps (Fig. 3 (a)) indicates that the 
system is initially in the oriented amorphous state. We 
readily notice reflections on the equator at q± — 14 nm^^ 
which are attributed to the nearest neighbor interchain 
length, and those on the meridian at gn — 52 nm^^ which 
correspond to half the fiber period of the polymethylene 
chain. At 4 ns. Fig. 3(b), broad layer streaks are visible 
along the layer line B,t q\\ — 26 nm~^, and the reflections 
at q± = 29 nm^^ on the equator are slightly sharper and 
stronger than in Fig. 3(a). It is clearly found in the 
snapshots that the degree of the order at 4 ns is consid- 
erably higher than that at 10 ps. At 8 ns. Fig. 3 (c), the 
layer streaks at g|| = 26 nm~^ become much more evi- 
dent, and the equatorial reflections around q± = 29 nm^^ 
come to sprit. Such diffraction pattern is quite similar to 
that observed in the hexagonal phase of PE, which ap- 
pears under high pressure and temperature p4[ or under 
high extension The snapshot in Fig. 3 (c) shows 

that the nearby chains tend to orient along the draw axis 
and form crystalline domains, which emerged around 4 
ns and gradually became lager. After 10 ns, marked in- 
tensity contrast appears on the streak at q^^ — 26 nm^^, 
and finally at 30 ns the diffraction come to show ordered 



crystalline pattern (Fig. 3 (d)) showing the development 
of high longitudinal order in the crystal. 

Around 15 ns the well-developed crystalline layer 
nearly filled the system. With further development of or- 
der, the amount of amorphous segments decreases. But, 
due to the periodic boundary condition, several parts of 
the layer surface contacted with the opposite layer sur- 
face of itself, which prevent the shrinkage of the MD cell 
along the draw direction and gives rise to the density 
decrease in the amorphous phase (Figs. 2 (b), 3 (d)). 

In summary, we have performed the MD simulations 
of the polymer crystallization from the oriented amor- 
phous states. We have succeeded, for the first time, in 
reproducing the fiber structure of stacked lamellae with 
ordered crystalline state. We have observed that, after 
initial relaxation for about 4 ns, the hexagonal order first 
appears and grows until around 10 ns, and then develops 
the highly crystalline state through longitudinal adjust- 
ments of the parallel chains in the crystal. Analysis of the 
fiber structure obtained, such as the conformation of tie 
chains and the structure of the interlamellar amorphous 
phase, and the detailed molecular process of fiber forma- 
tion are also very important problems, and they will be 
published in a separated paper. 
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